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Abstract 

In the paper we consider a hydrodynamical description of a nonideal Bose gas in 
one-loop approximation. We calculate an effective action which consists of mean field 
contributions and first quantum correction. This provides the equations of motion for 
the density and velocity of the gas where both mean field contributions and fluctua- 
tions are presented. To fulfill the calculation we make use of the formalism of functional 
integrals to map the problem to a problem of quantum gravity to benefit from a well- 
developed technique in this field. This effective action provides all correlation functions 
for the system and is a basis for a consideration of dynamics of the gas. Response func- 
tions are briefly discussed. Applications to the trapped bosons are reviewed. Together 
with Ref . 1^] the paper provide complete description of the condensate fraction and the 
deplition in the case of Bose condensed gases. 

1 Introduction 

Recently a lot of attention was attacted to the theory of a nonuniform Bose gas and its 
collective excitations. The interest was stimulated by the success of the experimental obser- 
avation of Bose-Einstein condensation systems of spin polarized magnetically trapped atomic 
gases at ultra-low temperatures |1|, 0, ^ and experimental studies of their collective proper- 
ties [0. On the other hand, new exciting problems such as a description of the evolution of 
Bose condensate from relaxed trap, dynamics of a collapse of the condensate for Li^ atoms, 
heating-cooling phenomena, various coherence effects for the condensate and so on are being 
posed both theoretically and experimentally. 
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Theoretically the evolution of the condensate and all related phenomena were investigated 
in number of papers H, ^ |l^ but all these considerations were formulated in Hartree- 
Fock (or mean field) approximation. It means that the calculations were based on Nonlinear 
Schrodinger equation (NSE) or Ginzburg-Gross-Pitaevskii equation (GPE) P]. It is obvious 
that the GPE describes strongly interacting Bose condensed gas well since it is possible to 
ignore non-condensed fraction. However, as it was noted in [P], this approximation fails 
to describe heating-cooling effects, fast dynamics or fast external potential variations, fluc- 
tuation effects around the edges of a trap and so on. To include the phenomena in the 
self-consistent consideration quantum corrections to the GPE should be calculated. 

Let us note that in the previous papers GPE plays double role. First of all it describes 
the Bose condensate fraction. But it also gives the dynamics of full density since in this 
approximation the deplition was ignored. Based on GPE collective excitations and the 



deplition were considered in RefsM, 11, 12, 12, 14| but the it was not self-consistent and the 



interaction of the condensate fraction with the noncondensed one have not been taken into 



account [|I6|. As a result, to leave the mean field approximation the GPE equation can be 
generalized in two different ways - we can write an equation for the condensate fraction as 
it was done in |P and its analogue for the full gas density. In the paper first (one-loop) 
quantum correction in the hydrodynamical approximation was calculated to generalize the 
GPE for the condensate fraction including effects of the interaction of the condensate with 
collective excitations (noncondensed fraction). To calculate the deplition in the framework 
we have to find Green functions for the corresponding effective action. Being principally 
straightforward this is not an easy technical problem. However there is a simpler way to find 
the diplition as a difference between full density and the condensate density if the first is 
calculated in the same approximation (it is not difficult to see that the mean field equation 
for the density is once again GPE as it was noted above). In the paper we solve the problem 
and give dynamical equations for the full density of the particles which generalizes GPE 
equation with an account of first quantum correction to the mean field picture. This allows 
us to calculate the deplition as well. 

The paper is organized as follows. In the next section we introduce functional integration 
to calculate the effective hydrodynamic action for the system and define one-loop quantum 
corrections which are responsible for the fluctuation effects. These quantum corrections are 
evaluated using (^-function regularization for functional determinants and an explicit and 
very simple expression for the one-loop quantum corrections is presented. This permits in 
Section 3 the derivation of the quantum corrections to the equations of motion for the gas 
in closed and explicit form readily used in numerical calculations. Here corrections to the 
Nonlinear Schrodinger equation (due to the interaction of the mean field background with 
the fluctuations) are obtained. Section 4 is devoted to a consideration of Green functions 
and response functions. It is shown that all Green functions in the 1-loop approximation can 
be obtained from the effective action of Section 3 by functional differentiation. We conclude 
the paper with a discussion of the possible applications of our results. The analysis does not 
depend on the details of the confining potential and can be used in a variety of problems. 
The only simplification is common hydrodynamic approximation. 
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2 Hydrodynamical picture 



As it was said above we consider nonideal Bose gas with point-like interaction in external 
potential. The Hamiltonian of the system then can be written as 

2m 



H 



dV 



where v{x,t) is an external potential. For the case of experimentally trapped atoms the 
potential is supposed to be harmonic and can be varied or even switched off at some moment, 
say t = 0. More precisely it means that for the case of magnetically trapped atoms v = 
0{—t)a]_^{t){p^ + ^'^{t)z^), = + y^, (a_L is oscillator lengths) has to be put in Eqn(|l]). 
Another parameter in Eqn(|l|) is I — the s-wave scattering length for atoms in the system 
18[| . The operators are Bose creation and annihilation operators. This Hamiltonian 

is a basis for the calculation of quantum mechanical and thermodynamical quantities. For 
example, the corresponding vacuum-vacuum transition amplitude is given by the following 
functional integral: 

= /D^+(x,r)DV^(x,r) exp[iS/h] , (2) 



with the action 



S 



dr 



'2,771 

Vip^Vip + {v - -^fi)i!'^ilj + AttIiP^iP^iIjiI) 



(3) 



where the parameter fi is the chemical potential, controlling the number of particles in the 
system. 

Now we will derive the hydrodynamic description directly from the functional integral 
formalism, which simplifies the description. We will consider zero temperature case but the 
generalisation for the case of finite temperatures is straightforward (though cumbersome) 
using Keldysh technique |ll9|. As it is explained in Appendix A, we are looking for the 
effective action which describes all physical quantities for the system. For example, the 
effective action provides all Green's functions in the same approximation used to calculate 
the effective action itself. We will obtain it in a hydrodynamic one-loop approximation. To 
clarify the description hydrodynamical variables, density and velocity, should be used. More 
precisely, we start with the action (pi) 



S 

h 



dr dV 



ihip^^^ 

OT 2m 



2777, 

h 



(4) 



and rescale variables as r ^ T^''"' ^ ^ following action form: 

_d'ijj 



S 

h 



dr dVi itp-^ 



dr 



V^+VV' + {v- + Ki^^ipf 



(5) 



(We will measure the action in terms of h everywhere below). Now we change field variables 
to the hydrodynamic ones: 



^{x,t) = Jp{x,T)e 



-i(p{x,T) 



ip+{x,T) = Jp{x,T)e 



J(P{x,t) 
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Then the action (|^) takes the form (up to a complete derivative term): 



(6) 



It is not difficult to see that the classical equations of motion for the action lead to the 
equations: 

^ - (V^)2 -v + ^-2lp + VP = , (7) 

+ 2V(V<^ ■ p) = , (8) 

which are equivalent to the Gross-Pitaevskii equation [|12] after the introduction of the 
velocity variable c = —2Vip instead of (p. In velocity-density variables these equations 
look as hydrodynamic equations for an irrotational compressible fluid: 

— + v{^- + 2v-2fi + 4lp-—V'^j = 0, (9) 

|^ + V(c-p) = 0. (10) 

Now let us shift our variables by the zero-order mean field solution (see Appendix A for 
details): 

p p + (7 , ip ^ (p + a 

and we will keep up only terms including the square of a and a terms. Then the action 
transforms to S{p, ip) + Si{p, ip, a, a) where the action Si has the form: 

/CO r I f)fy 

drJdV I —a - p{Vaf - 2aV ■ ^pVa - la^ 

"Ap 3(Vp)2 



p 2p2 



(111 



Last two terms vanish in the hydrodynamical limit since p is a large variable. Hence the 
action 5*1 takes form: 

^1 = p <\t j <W jl^fx - piWaf - la^ - 2aVip ■ Vaj . (12) 

Integrating in functional integral over the field a we get the action for a field only: 

The corresponding determinant is an inessential constant. 

Now we will use the large parameter p (for harmonic trap p ~ N'^/^l'^/^ if co-ordinates 
and / are measured in harmonic length units a±). Let us introduce the following large number 
Po = max IpI such that p = — ^ 1, t = Ty/AlpQ, v = ~ 1. Then 



dV A'^'d^ad^a. 



(14) 



where the matrix A has a form 



/ -1 



^^1 



V2 



A 



V3 \ 



Vi p-V{ -V1V2 -V1V3 
V2 -V1V2 P-Vl -V2V3 
\ Vs -V1V3 -V2V3 p-vl J 



with its determinant det74 = — p^. 

Our next step is to cast the action in the covariant form. To do this we introduce auxihary 
metric g^i, such that 
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y^det (11^^-11) 

One can easy to find the covariant metric from the equation above: 
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V-det(||^^1 



that leads to the form: 



/ -1 



-3/2 



Vl p-vl 



V2 
\ V3 



V2 
-V1V2 

-V1V2 p-vl 



V3 \ 
-V1V3 

-V2V3 



-V1V3 



-V2V3 p-vl I 



Finally we get the metric ||5'/ii/i|: 



\9^iv I 



^1/2 



^ -p + V^ Vx V2 V3 \ 

Vl 10 

V2 10 

V ^3 1 / 



with the determinant g = det = —p^- 

In this metric the action takes a covariant form 



dx J—gaOa 



dx 



ga- 



'-9 



^d^~g^^J-~gd,a . 



This is the central equation of the paper. 



Now the effective action is Ti = -tr ln[(po/4/)^/^n]/2 with □ 



(15) 



-gdy. Let 



us again make use the existence of large parameter in the system. Indeed, the multiplier 
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gives a possibility to express main contributions to the determinant such as (see 
Appendix B): 



tr ln[(po/4/)^/^n] = tr ln[n] 



-tr \n{po/Al){%{^) - L{a] 



~-tr ln(po/4/)($o(n) 

since in our regularization tr ln[n] is order of Seeley coefficient ^oC^) = /da;dr ^Z—g^^o. 
Returning to the initial variables all curvature tensors and the metric are written in the 
"physical" variables (i.e. without tildes). Moreover ln(po/4/) 3> ln(p) so that we can sub- 
stitute ln(p/4/) instead of ln(po/4/). Summarizing, we obtain an expression for the first 
quantum correction to the effective action: 



Ti = r; + r;' 



(16) 



where 



r; = -1 1 dxdr v/=^ln(p/4/)vl>o(n) 



with the metric 



r;' = - ln(po/4/)L(n) . 



vi 10 
V2 10 

V ^3 001/ 



and Vi = 2diip. 



3 Quantum corrections to equations of motion 

In this section we derive from (^) required quantum corrections to equations of motion (see 
Appendix A). To make the consideration self-consistent we should not differentiate ln(p/4Z) 
because in our approximation it behaves like a constant. In that case, the action we have to 
vary is covariant again, we find: 

dp dp bg^y 4 bp 



bip bif bg^y 4 b^p 
Now, 

= -iln(p/4/)-^ Ux ^g^o{^) = ln(p/4/)-A$„(n) 

IfMu 4 bgf,^ J 4 bg^y 
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Let us remember that 



$o(n) 



(47r) 



— RiivR^'^'^ 
180 ^ 



We note that one of the terms above is irrelevant: 

1 



j dx y/^V'^R = /dx ^/^-^d^^/^g'^'d^R = Jdx d^.i^^/^d^R) 

since this is just a complete derivative. Hence for our purposes it is sufficient to consider 
only 



<l>o(n) 



(4vr) 



da; 



-g 



72 



Rii,,R^'^ H R 

180 ^ 180 



J^pucrp 



Long but quite straightforward calculations lead us to the following result for the func- 
tional derivative (see Appendix C for details of the calculation): 
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5g 



$o(n) 



36(47r) 



-9i-ii 



-R^ 



-RapR"'' + -RapapR"'"^^ 

5 5 



+ RpuR — -R^iaRy" + -Rp.apaR'u"'"^ ^ Qfiu^R + 2^"^ u} R 



+ 



1 r 



+ -VPV'R^^p, 



taking into account 



5 



6 



-.crp 



$o(n) 



Sgpu Sg" 
This means that the equations of motions with 1-loop quantum correction take the form: 



V + n 



1 



5gpy 5T\ 1 



Jpv 



+ 4 ln(p/4/) 



1^ + 2V(V^ ■ p) 



= 0, 



0, 



(17) 



These two equations (|T7p and (|T^ replace the mean field hydrodynamic equations (|^J^) for 
the density and velocity of the Bose gas and are the main result of the paper. They contain 
contributions of fluctuations to the dynamics of the gas. For example, we note that the 
additional quantum pressure is due to dependence of the Bogoliubov particles determinant 
on the density of the condensate while local density changes (see continuity equation ([TsD ) 
comes from phase-dependence of the determinant. There is no difficulty to see that the latter 
correction obeys particle conservation law (that means it is a divergence) . 

To complete this section we consider the stationary situation when the equilibrium con- 
densate velocity is equal to zero there is no local particle transfer. However we think that 
equations (^) and (|T^ are of the main interest in non-equilibrium problems such as evo- 
lution of the gas under changing of trap shapes and so on. In equilibrium case of trapped 
atomic gases, the velocity field is equal to zero as are all time derivatives in the equations 
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of the previous section. This allows us to greatly simplify the equations and find quantum 
corrections in closed and compact form. Indeed, it is cumbersome but straightforward to 
check that in the regime there is no corrections to continuity equation (^), L{0) = 1, and 
there is only the contribution to the quantum pressure term in equation (P). This means 
that the equation of motion take the form: 

-j + v-fi + 2lp-—V^^+SP = 0, (19) 

where the correction to quantum pressure 6P has the form: 

6P = ( ^AHnp+^{d. In p) {d.A In p) 

48(47r)2^V 20 20^ '^^^ 

27 87 
+ -(9Alnp)(aAlnp) - -(Alnp)2 

6 T 

- — {di In p) {dk In p) {didk In p) 

1Q7 417 \ 

-(di In pfA In p - — - [{di In pf]^ (20) 



160' " ' 1280' 

where roman indices correspond to spatial derivatives. 

This implies the mean field approximation is applicable at the center of the trap where 
the correction 6P is negligible. On Fig.l we plot the contribution in the central region of the 
trap. Since the contribution is particularly important near to the classical turning points 
we do not make use the Thomas-Fermi approximation and exploit more correct analytic 
expression. The profile of trapped condensate was obtained in Ref. pT[ in the mean field 



approximation as an interpolation between small density and high density expansions; this 
demonstrates very good agreement with numerical solution of the corresponding mean field 



equation |22|. The expression for the profile with particles has the form: 



p(r) = i^iy [2/A^exp(4C - r^)] , (21) 

where W{x), defined as the principal branch (regular at the origin) solution to the Eq. 
We^ = X [||, is a standard MAPLE function H; Constant C in Eqn.(|l]) is to be deter- 
mined from the normalization condition / d^x p(r) = [|25[. The profile of the condensate 



density calculated from eq.(21) is shown in Fig.l for N = 10^ atoms. 

Two points should be noted. Firstly, the magnitude of the correction in the central 
region of the trap is of order of 10~^ relative to the main interaction term (see Figs 2) 
(the contribution in the intermediate region is shown on Fig. 3). Secondly, the sign of the 
correction term changes as can be seen from Fig. 4. This means that the character of exchange 
interaction due to excited states depends on the density of the condensate and leads to an 
effective attraction for density higher than some value and to an effective repulsion for a 
smaller density. However we have to remind that the leading contributions in the region 
are much larger than the corrections and after a self-consistent solution of the corrected 
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Nonlinear Schrodinger equation is performed the profile of the full density will not change 
significantly in this region of the trap. 

Since the sign of the contribution 5P became positive, the value of this term starts 
to increase and becomes quite considerable. However, at some point, our hydrodynamic 
approximation has to fail and we need other correction terms which will stop the increase. 
Although we cannot treat them properly it is possible to give an approximation for this 
region based on the low density expansion. Up to now we dealt with a high density expansion 
denoted dPh- It is easy to calculate the leading contribution in low density expansion that 
controls the behavior of depletion interaction correction at large distances. In section 2 we 
saw that this correction corresponds to sum of all the one-point diagrams for the system 
with the action @. In the low density regime the main contribution 5Pi is due to a one loop 
diagram. It means that the correction to eqn. (|T^ 5Pi is 



5Pi = Al5p . (22) 

It is hard to give exact expression for the deplition in this region since the comprehensive 
solution is required. We estimate the depletion as one fourth of the local density of the 
particles and use this estimate in Fig.5. It shows plots of SPi and SPh which let us estimate 
the contribution in the intermediate density region. It is easy to see that the two curves 
meet at point r ~ 6.13 where the density of the condensate p ~ 1 and both (high and low 
density) approximations fails. This however let us estimate the magnitude of the correction. 
It is obvious that the crack at the figure should be smoothed and this comes from high order 
corrections. We cannot calculate them explicitly and the only interpolation should be used. 
The magnitude of the correction to the Nonlinear Schrodinger equation is ~ 1 and in the 
region is compared with both kinetic and repulsion energies as it is shown in Fig. 6. 

We see that the only region where the corrections are not small is the boundary region. 
Here the corrections are order of both kinetic and interaction energies and can play significant 
role in the determination of the profile of the condensate close to classical turning points. 
The Thomas -Fermi approximation does not work in the region and hence both Stringari 
and Wu and Griffin WM approximations are inadequate. This explain the discrepancy 



of our results. 

Other point to note is that the form of Eqn.([T^) coincides with the form of the correspond- 
ing equation in the case of the condensate of the gas 0, i.e. there is only one generalization 
of GP equation for the stationary case. However, solutions for the condensate density and 
full density can be different since the generalized GP equations is highly nonlinear and its 
solution crucually depends on initial conditions. 



4 Green functions, response functions and the like 

Equations of motion and quantum corrections for them are the central questions of previous 
consideration. However in many applications related to experiment it is very important to 
analyze various response functions and form-factors. They can be expressed in terms of 
some combinations of Green functions of the theory. That is why in this section we shortly 
consider a calculation of Green functions in the effective action formalism. 
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As it is shown in Appendix A the effective action (in any self-consistent approximation) 
allows the evaluation of all Green functions in the same approximation by just taking of 
variational derivatives. Since the effective action in 1-loop approximation was obtained the 
problem of Green function calculation is the problem of functional differentiation only [^. 

We now formalize all said above and give formulas to evaluate Green functions in effective 
action formalism. Let W{J) being generating functional of connected Green functions. Then 
quantities 

».4,,,...,,„) = _i_..._J_»V) (23) 

are connected Green functions in the external field J. The effective action is defined then 
by the Legendre transformation of W: 

T{a) = W{J{a))-aJ{a) , « = (24) 

oJ 

where the function a{x) = {(p{x)) = Wi{x; J) is the first connected Green function (pSj). It 



is easy to see that the functions a and A are related by the second relation (^). Then the 
quantities 

s s 

define IPI (one-particle-irreducible) Green functions for the theory with mean field a{x) = 
{ip{x)). Knowledge of IPI Green function is equivalent to knowledge of any (whole or con- 
nected) Green functions for the corresponding system. For the first IPI Green function we 
have from (E^) 

ri(x) = -J(x) . (26) 

All connected Green functions can be expressed in terms of IPI ones. Indeed, differentiating 
(^) J we obtain 

W2T2 = -1 , W2 = -T^' (27) 

or, in expanded form, 

JdzW2{x,z)T2{z,y) = -5{x~y) . (28) 
Differentiating now the second relation in ( P7| ) on J and using the following rule 

5J 2 

we can derive expressions for all the higher connected Green functions through the IPI ones. 
For example, for the third connected function we have 



1 3 

2 



F 



3 



or 



W3{xi,X2,X3) = - J dyidyidys F2 ^(si, ?/i)F2 ^(x2, l/2)r2 ^(xs, y3)F3(?/i, 1/2, 2/3) (29) 
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and so on. It means that since we know IPI functions we have to solve the differential 
equation (pSD for the two-point correlation function and then find all other connected Green 
functions by integration of IPI functions with two-point correlators as in Eqn.(^). 

Putting J = in (^) we get equations of motion for p and ip ([T7| , p!8|) . To obtain, 
for example, the Green function {{p{x) — {p{x))){p{y) — (p(y)))) we twice differentiate the 
effective action on p, substitute solution of equations of motion and finally invert the result 
in the sense of the kernel of a integral operator. 

As an example of usage of this technique we calculate two-point Green function in the 
mean field approximation (tree or 0-loop approximation). In this approximation the effective 
action has the form: 



{v - p)p - Ip^ 



(30) 



such that, in the equilibrium hydrodynamic picture, the matrix of second variational deriva- 
tives of the effective action can be written as 

I -01 - — 

dt 



5'T 



6p{x)6p{y) 6p{x)6(p{y) 
6^T 6^T 



-21 
d 



\ dt 



2V-pV 



V 6^{x)6p{y) 6ip{x)6ip{y) J 

where p = (p). The matrix of two-point connected correlation functions of fluctuations 

p{x) = p{x) — {p{x)) and (^(x) = 0{x) — {(p{x)) 

(p{x)p{y)) (p{^)!p{y)) 

{^{x)p{y)) {Hx)Hy)) 

is defined solution of the equation 



G 



I . 



Solving the equation we get the following expressions for the correlators 

1 d 



{p{x)p{y)) = 6{x-y) 



21 dt 



{ip{x)p{y)) 



1 d 



{p{x)^{y)) = -^^(^(^)^(?/)) 
and the equation for the phase-phase correlator {(f{x)ip{y)): 

[ 21 dt^ 



2V ■ pV {^{x)^{y)) = 5{x - y) 



(here x denotes both space and time variables). It means that the phase-phase correlator is 
the Green function of the wave equation. 

In summary, various response functions, form-factors and Green functions are evaluated 
from effective action in any given approximation taking variational derivatives and integrat- 
ing them with the two-point correlation function. This can be done in general using formulas 
of this section. 
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5 Conclusion 



In the paper a self-consistent approach to the calculation of 1-loop quantum corrections to 
the Gross-Pitaevskii equation for a Bose gas density due to fluctuations around mean field 
configuration are calculated. To do this the hydrodynamic approximation was used. In this 
approximation excitations are equivalent to sound waves on a mean field background. This 
opens the possibility to use methods of quantum gravity and theory of quantum gauge fields 
where the problem of calculation of effective action for gravitational and gauge backgrounds 
(or more precise, quantum corrections due to other quantum fields) is common problem. 

Many of the methods to approach the problem were developed during the last four 
decades. They are (^-function regularization for determinants of operators, Schwinger-De- 
Witt-Seeley expansion of heat kernels ^7\, , covariant methods of calculation of Seeley's 

coefficients and the covariant perturbation technique for effective actions |^0|, Although 
the methods are well-known in ffeld theory they are not so familiar in condensed matter 
physics and the theory of coherent systems. An accurate account of corrections is very 
complicated even for a few ffrst orders and the solutions are obtained by using the covariant 
perturbation technique and curvature expansions. On the way effects nonlocal in space 
and time appeared. In quantum gravity framework such effects play significant role in the 
gravitational collapse problem. Hawking radiation and so on. 

However to benefit from it sometimes very general and only basic information about 
principle scales in systems in question is required. Indeed, this is the only information needed 
to extract leading logarithmic contributions while the calculation of other corrections is much 
more complicated. That is why we think that the developed approach can be widely used 
in many problems which have not much to do with Bose-condensation of trapped atoms or 
liquid Helium. 

At the end let us stop on other applications of the presented method. It is easy to imagine 
other condensed matter examples of problems treatable by the same technique. Clearest and 
simplest of them are antiferromagnets, Josephson arrays and superconductors in nonuniform 
(in space and time) external magnetic fields. However many of other physical systems in a 
non-equilibrium background are potential field of applications. 

Returning to Bose condensation, the paper together with our previous work |Q where 
the generalization of Gross-Pitaevskii equation for the condensate fraction was given allows 
to consider the condensate fraction and the deplition on the same basis since the deplition 
is a difference of the full density and the condensate density. However, the results of this 
paper can be used even without refering to the Bose condensation phenomena. 
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Appendices 



A Effective action 

In this Appendix we state some basic information concerning the formahsm of an effective 
action following Let us consider a set of fields = ip^{x) = {ip{x),ip~^{x)}. The 

classical dynamics of these fields is described by the classical equations of motion 

(A.1) 



derived from the classical action functional 

S{^)= [dx C{ip{x),dM^)) ■ (A-2) 



After quantization the fields become field operators (p acting on some state vectors of Fock 
space. Let us suppose that there exist some initial |in) and final (out| vacuum states defined 
in some appropriate way. The vacuum-vacuum transition amplitude can be presented as a 
Feynman path integral 

(out|in) = Z{J) = Idv? exp 1^ (S{^) + Ji^)^ (A.3) 

Here we introduced the classical sources Jj to investigate the linear reaction of the system 
on the external perturbation. All the content of quantum field theory with all quantum effects 
is contained in the following Green functions 

\ (out|T((^*i ■■■0'") |in) , 

where T means the chronological ordering operator, i.e. the fields must be arranged from 
left to right in order of decreasing time arguments. 

It is easy to show that all these Green functions can be obtained by the functional 
differentiation of the functional Z{J), that is called, therefore, the generating functional. 
Moreover, factorizing the unconnected contributions one comes to generating functional for 
connected Green functions W{J) = — ih In Z{ J) 

The lowest connected Green functions have special names: the mean field 

V') = = Jjr (A.6) 



and the propagator 

(((^' - 0^) - 0')) = -thG^\j) = -^h^^ (A.7) 
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Further, the connected Green functions are expressed in terms of vertex functions. The 
generating functional for vertex functions is defined then by the Legendre functional trans- 
form 

r((/.) = w{jm Jki^)^ (A.8) 

This is the most important object in quantum field theory. It contains all the information 
about the quantized fields. 

1. First, one can show that it satisfies the equation 

r.(0) = ^ = -U<t>) (A.9) 

and, therefore, gives the effective equations for the mean field. These equations replace 
the classical equations of motion and describe the effective dynamics of background 
fields taking into account all quantum corrections. Thus r(0) is called usually the 

effective action. 

2. Second, it determines the full or exact propagator of quantized fields 

-r,^,^^"^^," (A. 10) 

where 5f = 5^5{x — y), and the vertex functions 

r„ = r,,...,„, (n>3) (A.ll) 

This means that any iS-matrix amplitude, or any Green function, is expressed in terms 
of propagator and vertex functions that are determined by the effective action. 

3. At last, when the test sources vanish the effective action is just the vacuum amplitude 

(out|in)|_^^^ = exp(ir/^)|^^o (A.12) 

Let us rewrite the definition of the effective action. To define the path integral it is 
convenient to make a so-called Wick rotation, or Euclidization, i.e. one replaces the real time 
coordinate to the purely imaginary one t —it and singles out the imaginary factor also from 
the action S — > —iS and the effective action F —iV. Then the metric becomes Euclidean, 
i.e. positive definite, and the classical action in all good field theories becomes positive 
definite functional. So, the Euchdean effective action is defined to satisfy the equation 

exp [Im) = /dV' - - } (A.13) 

This path integral is still not well defined. There is not any reasonable method, except 
for lattice theories, to calculate this integral in general case. The only path integrals that can 
be well defined are the Gaussian ones. Then the full path integral can be well defined as an 
asymptotic series of Gaussian ones. This is just the quasiclassical, or WKB, approximation 
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in the usual quantum mechanics. We decompose the fields in the classical and quantum 
parts 

ip = ^+Vhh (A.14) 

and look for a solution of the equation for the effective action in form of an asymptotic series 
in powers of Planck constant. 

r(0) = s(0) + ^^"r(„)(0) (A.15) 

n>l 

Then all the coefficients of this expansion can be found. They are expressed in terms of the 
well-known Feynman diagrams. The number of loops in these diagrams correspond to the 
power of the Planck constant. We will be interested below in the so called one-loop effective 
action 

1(1) = — InDetF (A.16) 

where F^k = S^ik- 



B Determinants of elliptic operators 



Although Eqn.( [ATT^ ) seems to be very easy it is still ill defined. The point is it is divergent. 



This is just the well-known ultraviolet divergence of the quantum field theory. Indeed, we 
can rewrite the functional determinant as 

Indet A = Tr InA = ln[]A„ = ^lnA„ (B.l) 

n n 

where A„ are the eigenvalues of the operator A. This series is easy to show to be divergent. 

That means that one needs a regularization. This point was investigated very thoroughly 
by many authors and it is found that in quantum gravity and gauge theories the most 
appropriate regularizations are the analytical ones. The functional determinants can be well 
defined in terms of the so called (^-function. 

At first let us use a Wick rotation to produce a Laplace operator from the wave operator. 
It allows to make use of methods for the evaluation of the determinants of elliptic operators 
Now we are ready to introduce the C-function of an elliptic operator A: 



as,A) = Y.K' 

i 

where {A^} are eigenvalues of the operator A. Then 

TrlnA = --^as,A) 
To study the behavior of the ^-function it is common to use the formula 

oo 

Ar = ^^/dt f-'expi-tX,) . 
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Summing over all nonzero eigenvalues of A we get 

oo 

as, A) = — - |dt f^\Tr exp(-tA) - L{A)) (B.2) 
^ ^ 

where L{A) is a number of zero-modes of A. Representation (|B.2|) is valid when the integral 
converges. For nonnegative operator it always does converge as t — +00. Conditions of the 
convergency of the integral as t —* +0 depend on details of the operator A. For Laplace 
operator in 4D it converges as t — +0 if Re s > 2. 

There exists the well-known Seeley expansion for the Tr {exp{—tA)): 

Tr exp(-M) = ^ $_fc(A)r'= + p{t) , 

where is bounded by a constant times t as t ^ +0. $_A:(y4) are called Seeley coeffi- 

cients and play important role in the investigation of elliptic operators and their topological 
properties (for example in the Index theory). For example, for 4D Laplace operator = 
A; > 3, and 

$_2 = (47r)-2 , 

$_i = — (47r)~^|i? , where R — scalar curvature , 

and 

*» = (4.)- (-^V^fl + ifl^ - + • 

Here R^i^ap and Rfj,i, are Riemann and Ricci tensors correspondingly. They are defined by a 
metric of a curved space. The latter coefficient $0 is particular important for the calculation 
of determinants. 

Splitting (P3.2|) into an integral over [0, 1] and one over [1, +00] we get 



/OO rl 
dt Tr exp{-tA)t'-^ ^ J ^'^ pW^'"^ 



The singularity at s = turns out to be removable since lim^^o sT{s) = 1 and ({s) is defined 
by analytical continuation from the half-plane Re s > 0. 
From the relation Indet A = — C'(0) we find that 

IndetA = ^^^^ + r(l)($o(A)-L(A))- r^Trexp(-M)- 

A;>0 Jl t 

- J^'j(TTexp{-tA)-J2^-kiA)t-'^^ . 

This is ^-regularized In det A which we use in the paper. 

To conclude this section we give the relation between In det A and In det a A where a is 
a number parameter. It is easy to see that 

C{s,aA) = C{s,A)-a-' 
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because Xi{aA) = aXi{A). This leads to the relation: 

In det aA = In det A + In a ■ C(0, A) = In det A + In a • ($o(A) - ^(A)) (B.3) 
which we use intensively in the paper. 



C Variational derivative of the effective action 

In this appendix we give details of the calculation omitted in section 3 to do not interrupt 
the main stream of the consideration. To make the calculation more efficient we make 
use covariant form of the quantum correction since there exists well-know way to simplify 



covariant variation calculation p4 |. 

We are interested in the following variational derivative 



6 



Sg 



fii/ 



where 



$0(n) = — L- / dx ( — R^ - —RuuR"" + —RauapR^""""^ 

ov ; (47r)2i ^ ^V72 180 ^ 180 ^ J 



(C.l) 
(C.2) 



For variation of the volume element we have 

1 



5^g 



where A'^'^ is defined by as 



-9 



-Sg.ug'"' 



-9t^uSg 



111/ 



A' 



9 ■ 9 



and the following relation was used 



-9tiuSg 



fii/ 



For variations R^, R^^R'^'^, R^y^pR^"'^'' one can obtain 

5R^ = 26R-R = 26g^^RuR + 2g^^6Ri2 ■ R 
6Ri2R'^ = S{g'^g^^RuR34) = 26g^^Ri3R2^ + 26RuR^^ 
6R,2mR^^^^ = 5{g^''g^^g'-'g^^R,2Z4R^,7s)=A5g^^Rr^,^R-^''^'' + 25R,2uR 

where numbers 1, 2, 3, . . . mean indices fii, fi2, f^s, ■ ■ ■ respectively. Hence, 

1 



1234 



36(47r) 



dx 



-9 ^9- 



flU 



23f^^ \ ~ ^-Rl2-R^^ + --Rl234-R^^^^ 



2 4 

+Rjj,vR — -^RfiiRk + '^Rfj.usRi, 



123 



36(4 



g^^5Ri2R - l5Ri2R^^ + ^5Ri2^^R^^^^ 
5 
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So we need to calculate last three terms. 



For variations of Ricci and Riemann tensors we have [B4 
6R 



1 



6R 



-gf,i6g^'^R2uap + g^^lVp{6Tl^) - gf,iVa{STl ) 



1 r. 



Let us consider covariant divergence of some vector T'^: 



/xp I V (7 V p,^yup 
k 



(C.3) 



(C.4) 



'-g 



Hence 



is a pure divergence and can be dropped under integration. Using this we can easily derive 
the following expression for variations of the integrals we are interested in 

1 



J dx ^g^^5R^,R = Jdx V^Sg^" [-^{'^^., V.}i? + g^.OR 



dx ^^dRauapR^''''^ = / da; V^^g 



Rpi23K''' + 2V'V'R^ 



Summing all terms we obtain 



6 



6g 



pu 



36(47r) 



-g\ -i^gt^v 



pl2v 



-K — -KcrpK + -KcrpaiS-n. 



2 2 1 

+ RpuR — -RpcjR'u" + -^RpapaR'/''"' — gpu^R + 2^'^^^' '^ujR 



1 r 



+ 5 



-2V'^V.i?p. + BR^, + gpuV'^VPR.p] + ^V^V'^R^.pu 
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Figure caption. 



1. Fig.l Bose condensate density profile obtained in mean field approximation. The curve 
is drawn using the analytic Ansatz of Ref.[^. Lengths are measured in units of the 
oscillator length a_i_. The multiplier is defined as A = 47cl/aj_. 

2. Fig. 2 The calculated one-loop quantum correction to the Nonlinear Schrodinegr equa- 
tion in the central region of the isotropic trap. Lengths are measured in units of the 
oscillator length. is the number of particles (for the picture = 10^) and 6P is 
defined by eq. (pOD . 



3. Fig. 3 The calculated one-loop quantum correction to the Nonlinear Schrodinegr equa- 
tion at the intermediate region of the isotropic trap. Lengths are measured in units of 
the oscillator length. A^ is the number of particles (for the picture A^ = 10^) and 6P is 
defined by eq. (^ . 



4. Fig. 4 Change in sign of the one-loop contribution in intermediate region of the trap. 
Lengths are measured in units of the oscillator length. A^ is the number of particles 
(for the picture A^ = 10^) and 6P is defined by eq.(pO|). 

5. Fig. 5 Merging of High density approximation for the quantum correction (solid line) 
and the corresponding Low density one (dashed line). Lengths are measured in units 
of the oscillator length. In the region to the left to the point r ~ 6.13 High density 



approximation produces the correction (|20|) . In the region to the right to the point 



r ~ 6.13 Low density approximation ([2^ ) for the correction should be used. In the 
region of the boundary point both approximations fail and other quantum contributions 
should be kept. 

6. Fig. 6 Comparison of various terms in the equation for the condensate density. The in- 
terpolating quantum correction (solid line) of Fig.5, nonlinear interaction term (dashed 
line) and kinetic energy (circles) are pictured. Lengths are measured in units of the 
oscillator length. 
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Fig 



.3. Quantum correction: intermediate region 




Fig. 4. Quantum correction: change of the sign 

tI n5p/100 



Fig. 5. High and Low density expansions 




Fig . 6 . 



Comparison of contributions 




